Analytical form of light-ray tracing in invisibility cloaks 
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Abstract 

In this paper, we review the methodology of transformation optics, which can construct invisi- 
bility cloak through the transformation of coordinates based on the form invariance of Maxwell's 
equations. Three different ways to define the components of electromagnetic fields are compared 
for removing some ambiguities. The analytical expressions of light-ray and wave-normal ray are 
derived in spherical and cylindrical ideal invisibility cloaks created with any continuous radial 
transformation functions, and their physical interpretation is also given. Using the duality prin- 
ciple in anisotropic media, we prove that light-ray vector satisfies "ray-vector eikonal equation" 
corresponding to the usual "wave- vector eikonal equation". The results interpret why the wave 
vector maps to the ray vector transferring from the virtual space to the physical space, but not the 
wave vector. As an application, we investigate the special transformation functions which make 
the light-ray function satisfy harmonic equation. 
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I. INTRODUCTION 



Recently, invisibility cloak has attracted much attention and widely study its 
theory is based on the form invariance of 3-dimensional Maxwell's equations under a coor- 
dinate transformation [l|, 0-0] • Briefly, the property of material, such as e, fj, in a vacuum 
Euclidean space with an original orthogonal coordinate, will be transformed into a different 
form through a coordinate transformation. However, the form of Maxwell's equations is 
invariant. The transformed s and \x could be alternatively interpreted as a new different set 
of material properties in the original orthogonal coordinate system. Based on this point, the 
developments have generated a new subject of transformation optics Meanwhile, a con- 

formal mapping method has been applied to design a medium which can achieve invisibility 



in the geometric limit [2[. The validation of the invisible cloak designed with transformation 
optics method has been verified through several ways, including ray tracing approach [4, 17], 



full- wave simulations [lOj . and analysis based on Mie scattering [11-16]. The early studies 
in the invisibility cloak were mostly based on linear transformation proposed by Pendry et 
al. llMl3l]. Then different transformations have also been discussed [14|, ll6|]. Luo 



et al. proved that the permittivity and permeability created with a most general form of 
transformation can realize spherical invisibility cloaks and give the analytical calculation of 
electromagnetic fields for any spherical cloak using Mie scattering method [l5| . 

In the pioneering work of Pendry et al.[l], it had been noted that the Poynting vectors 
acting as the tangent vectors of light should propagate to surround the hidden area in the 
cloak layer. They used the canonical equations of a double anisotropic material optical 
system to describe the light propagation 4J]. According to the analogy between the curved 
space and double anisotropic media, Ref. [17J described the ray traces with the geodesic 
equation in the virtual curved system, and gave the proof that it is equivalent to canonical 
equations. Since the original Euclidean space (virtual space) is a vacuum, the Poynting 
vectors and wave vectors of a beam of light traveling in it are parallel to each other. Whereas, 
when mapping to physical space where a cloak is placed, Poynting vectors and wave vectors 
are not paralleling any more, which has been pointed out and discussed by Leonhardt in 



20l | . In this paper, however, we will expound this point from a different perspective by using 



"ray-vector eikonal equation" and corresponding relations. 

The paper is organized as follows. At the beginning, we compare three common ways 
to define the components of electromagnetic fields in 3-D curvilinear coordinates, which are 
named Minkowski's definition, Landau's definition, and the components in normalized bases. 
To do this is benefit to clarify some ambiguities caused by the mistiness of definition. After 
that, We repeat the fundamentals of transformation optics, which contain the corresponding 
relations between the curvilinear coordinate system of virtual space (S system) and real 
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physical space (P space), and emphasize that the validity of transformation optics is not 
only based on the form of Maxwell's equations, but also based on the correspondence of 
boundary conditions between the two systems. 

After explication of the fundamentals, the analytical solution of electromagnetic fields is 
derived in spherical and cylindrical cloaks with arbitrary radial transformation functions. 
The Poynting vectors and wave vectors are further calculated, and the analytical expressions 
of light rays and wave-normal rays are obtained, respectively. Then, we apply the expressions 
to four types of cloak made by four different transformation functions. Furthermore, we 
provide the physical interpretation to the analytical solution of light-rays, and verify it 
satisfies the geodesic equation in S system discussed in [17] ■ Through duality principle of 
anisotropic material, we obtain "ray-vector eikonal equation" corresponding to the usual 
"wave- vector eikonal equation". And through the corresponding relations, we prove that 
the covariant component of wave vector hi in S system maps to wave vector /c (P)i in P space, 
however, its contravariant component k % maps to ray vector s* p) in P space. This is the 
precise reason that wave vector and Poynting vector are split in P space. 

Finally, as an application, we find out the general form of transformation functions that 
makes the light-ray expression satisfy harmonic equation. It is just the conformal transfor- 
mation function, which Leonhardt applied to design invisibility device in the geometric limit 
In this case, Poynting vectors and wave vectors have neither divergence nor curl. 



II. COORDINATE TRANSFORMATION OF ELECTROMAGNETIC QUANTI- 
TIES 



A. Definition of electromagnetic fields in 3-D orthogonal curvilinear coordinates 

For an arbitrary 3-D curvilinear coordinate system {x 1 } to describe spatial part of a 
lat spacetime, its bases e*j = dr/dx % are often called coordinate bases or holonomic bases 
2l| . and its spatial metric is 7^ = e*j ■ e*,-. The components of a tensor T are defined by 
f = T^eiCj. 

The component forms of Maxwell's equations in any 3-D curvilinear coordinate system 
are: 

1 f) R' f) R* 1 

VtB* = —dii^yB') =0, — + ji%E h = — + —e ijk djE k = 0, (la) 

y/l dt dt yfi 

1 _ (jjji _ 8D i 1 

Vi& = —d^D 1 ) =p, -— + e i i k d j H k = -— + —e ijk djH k = 0, (lb) 



where Vj is 3-D covariant derivative, 7 is the determinate of 7^, e y is 3-order completely 
antisymmetric symbol, and e yfe = e* jfe / ^7 is 3-order Levi-Civita tensor. We use Latin indices 
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i,j, • • • to denote 1 to 3 for three spatial components. The components of electromagnetic 
vectors in these equations are resolved in holonomic bases, and this mode of components is 
named Landau's definition (see Appendix A for general Landau's definition in 4-D curved 
spacetime). We add a bar on these components to discriminate Minkowski's definition, 
which will be introduced below. 

In this representation, E, H, D, B are all treated as 3-D spatial vectors, which imply 
their components satisfy the law of vector transformation T l = (dx l /dx' k ')T' k \ meanwhile 
e and jl are treated as 3-D spatial tensors, and their components satisfy the law of transfor- 
mation e lj = A],A 3 j, i' 1 ' 3 ' , fi lj = A^Ay , where A\, = dx l /dx' 1 ' is the element of Jacobian 
matrix. In vacuum, e^'/eo = ft 1 -* / fJ-o — Y J ■ 

In transformation optics, a different way to define electromagnetic fields in curved coor- 
dinate is established in order to hold the expressions of Maxwell's equations in Cartesian 
coordinate . If we define electromagnetic quantities as 



D' 



B' 



(2a) 



e« = yfy? j , = y/lP?- (2b) 

P = f = (2c) 

substituting them into Eqs. (PQ), the form of Maxwell's equations in curved coordinates 
changes to the form appearing in Cartesian coordinates. As a result, the form of Maxwell's 
equations expressed by ordinary derivative is invariant from coordinate transformation. 
We name this mode of definition as Minkowski's definition (see Appendix A for general 
Minkowski's definitions). 

Note that the above components D l , B l , e % \ /i*- 7 are not the real components of the 
corresponding vectors or tensors in the curvilinear coordinate, but they can be regarded as 
the components of y/jD, y/jB, y/jjl respectively. These quantities are all vector or 

tensor density, which satisfy the law of transformation 



D' 



det (A\,y l 



A), D' 



B' 



det (A*,)- 1 



A;, B' 



det (A*,) 



-i 



a;,.v;, e' ij 



det(Ai,) _1 A*,AJ,^. 



In vacuum, the permittivity and permeability have the simple forms as 



r'J 



'j 



(3a) 



(3b) 



(4) 



The components resolved by unit basis are also commonly used to expressing tensors in 
curved coordinates. The set of normalized bases is defined as ii = ei/\\ei\\ = ei/ y/ju, yet 
these bases are anholonomic in that we usually can not find a set of coordinates to satisfy ii = 
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dr/dx 1 . The components in orthonormal bases are also defined by T = e^ej, and there 
is no distinction between covariant and contravariant components when the coordinates are 
orthogonal. For a 3-D orthogonal curvilinear coordinate system the set of orthonormal bases 
is corresponding to the tetrad in 4-D spacetime (see Appendix A). Since all the components 
of a vector in unit bases have the same dimension, they act as the measurement of the vector 
in each direction. In terms of the relation of holonomic and anholonomic bases, we can get 

e {ll) = hihj^ = b*bL e v ^) = hihj ji» = (5) 

where hi = y/ju- For orthogonal coordinates, ^7 = hi h 2 h s . In Eqs. (JSJ), we have canceled 
the Einstein summation convention. For vacuum, s^/sq = fi^/fio = 5ij. 



B. Three coordinate systems in transformation optics 

The method of transformation optics starts from a virtual 3-D flat Euclidean space with 
medium of vacuum. An orthogonal coordinate system (Cartesian or curvilinear coordinate), 
S' system, is employed to describe the virtual space. In S', the coordinates are written 
as x' 1 and the metric is 7'j/j'- However, another coordinate system S with coordinates 
x l and metric 7 ij can be established to describe it. The transformation of coordinates is 
x l = x l {x' % '). 

If we investigate electromagnetic fields propagating in the virtual space, Eqs. (BJ give the 
Minkowski's form of permittivity and permeability describing vacuum either in S or in S'. 
And Eqs. ([3]) present the transformation law of Minkowski's form between S and S' system. 

Due to the invariance of Maxwell's equations under Minkowski's definition, we can alter- 
natively interpret the Minkowski's form of permittivity and permeability in S system as the 
properties of a real material in physical space (P space). The P space is also an Euclidean 
space and described by the same orthogonal coordinates as S'. We use x (P) l to denote the 
coordinates of P, and also add the subscript "(P)" to other physical quantities in P space 
for definitude. Since we employ the same coordinate system as S' system, the spatial metric 
7 (P)ij . of P space hold an identical form with j'^/jt- On the other hand, the Minkowski's 
form of electromagnetic quantities in P space takes the same formula as that in S system 
of virtual space (see Fig. 1). 

In previous articles, the S system of virtual space is often called a curved space because 
of its nontrivial metric. In mathematics, it is certainly a metric space or manifold. However, 
considering its zero curvature, we only treat it as a curvilinear coordinate system of virtual 
flat space. In this paper, we only call a system space when they have a real physical 
background of spacetime, such as virtual space and P space. But for the coordinate systems 
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FIG. 1. Relations among Physical space, S' coordinate system and S coordinate system of virtual 
space. 

used to describe a space, we only call them system but not space. Nevertheless, we do not 
distinguish the P space itself and its coordinate system, because we only set one system to 
describe it. 

The corresponding relations of electromagnetic quantities between physical space and S 
system are shown in Table HI There are several points to emphasize. First, corresponding 
means the Minkowski's components of these quantities regarded as functions of coordinates 
have same mathematical forms in the two systems. Second, for an individual quantity, only 
the covariant (or contravariant) components hold the same form in the two systems, yet the 
other kind of components isn't congruent. For instance, DlpJx l , p S) and D l (x % ) have the same 
mathematical form, but the forms of D^ P)i {x\ p C) and Di(x l ) are different. Third, the identity 
is only valid for the components defined in Minkowski's method, however the components in 
Landau's definition and in the unit bases often take the different form in the two systems, 
unless the three kinds of components reduce to equivalent under Cartesian coordinate of 
physical space. 

TABLE I. Corresponding relations. 

Electromagnetic quantities 
Physical space : D l (p) E {P)i B\ p) H {P)i e (P) l - 
S system : D i Ei B l Hi e ij 



1 1 1 



Now our discussion is restricted on the condition that S system is an orthogonal coordi- 
nate system with 7^ = h^ij. Based on the corresponding relations, Eqs. (jSJ), and flSJ), we 
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derive the permittivity and permeability in physical space: 

H-y /i 2 /^(p)i h( P )j 



hp){ij)/ £ o — Vmiij) / 



E 



■f h (P) i h {P) 2 /^(p)3 (^y) z 



det (A*,)' 



MM. 



(6) 



Here we take the transformation matrix as a function of x l and x' l \ i.e. K\, = h\,{x % , x n '), 
meanwhile the independent variables should be x l {p) and x H in the above formula. 

For spherical invisibility cloak, if we select S' system as spherical coordinate system, 
define the inverse transformation from S' to S system as r' = f(r), 8' = 8, = 0, and 
insert it into Eq. ([6]), the constitutive parameters become 15] 



/ f 2 (r P ) \ 

(e iP){ij )/e ) = (n (P)(ij) /vo) = diag (P) f (r (P) ), /'(r (P) ) . 

V r (p) J v(p)) J 



(7) 



For arbitrary transformation function, it can be completely invisible, as long as f(a) = 
and f(b) = 6,where a, b are the inner and outer radii of the cloak layer, even if /(r) is 



15]. Similarly, the constitutive parameters for cylindrical 



continuous but non-differentiable 
cloak obtained by the mere radial transformation r' = f{r), 8' = 8, z' = z are 16] 

f( r (P)) W(?>)) /( r (P))f ( r (P))x 



(s iP){ij) /s ) = (fi (P){ij) /fi ) = diag 



W( r (p))' /( r (p)) 



(p) 



where the components are described by cylindrical coordinates. The arbitrary /(r) will 
achieve completely invisible only if it satisfies the boundary condition f(a) = and f(b) = b 
( see the proof in Appendix B). 



C. Boundary conditions for perfect invisibility cloak 

The congruity of Maxwell's equations with the corresponding relations between S system 
and P space is regarded as the foundation of transformation optics. However, despite of the 
existence on the same form of Maxwell's equations, they won't give the same solution unless 
the boundary conditions of the two systems are also identical. 

For example, the realization of spherical and cylindrical invisibility cloak requires the 
inverse transformation function satisfies f(a) = and f(b) = b. The f(a) = causes the 
domain r' < a in S' system to shrink into its original point in S, so that the coordinates 
with radial component r < a would never appear in the expression of electromagnetic 
fields in S system. It leads to decoupling of fields between the hidden region r (P) < a and 
cloak layer in the P space. Pendry et al. pointed out that the condition f(b) = b makes the 
constitutive parameters satisfy perfectly matched layer (PML) conditions £( P) (ed) = £ (p)(<^> = 
l/£(p)(rr) and AV)(6>6») = /V)(#> = V/VK"-} at t ne outer surface of the cloak (r (P) = 
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therefore the cloak is reflectiveness. However there is a more profound reason produced by 
the corresponding relation of the boundary conditions between S' system and P space. 

First, we consider a spherical cloak designed by general transformation x l (r f , 9', <$') from 
S' to S system. In P space, the boundary conditions at the interface r = b are 

(E (P) \ r=b+ -E m \ r=b _) x e r = 0, (D(r)\ r=b+ ~ D (P) \ r=b _) ■ e r = 0, (9) 

where e r acts as the normal vector of interface. The boundary conditions indicate the 
tangential component of E and normal component of D are continuous across the interface, 
and the two boundary conditions can completely determine E (P) \ r=b -, if E {P) \ r=b + is given. 
Since the metric of P space is continuous, the equality of components is always tenable no 
matter what definitions. Thus, we have the continuous relations in Minkowski's form 

E( P) e\ r=b + = E {P) e\ r=b _, £W| r=6+ = -E (P)1? i| r=6 _, D r (p) \ r=b+ = D r (p) \ r=b _. (10) 

Non-reflection means the electric field outside the cloak only has the incident wave. Ac- 
cording to the corresponding relations in Table HJ the electric field in S system should also 
satisfy the relations 

K\ r/=b + = E e\ r=b + = E e\ r=b - = A e E i'\ x 'i'( b -,e^y ( lla ) 



^ \r'=b+ \r=b+ 



E e 




= E e \ 




r=b+ 


" 1 r- 


Ecj, 




= E Ar-- 




r=b+ 




+ = 


-- D r 


r=b~ 



K'\r'=b+ ~ E Ar=b+ ~ E Ar=b- ~ A * E '%' L«' (6" ,6,0) ' ( llb ) 

det(Ai,)" 1 (He) 



where E' e ,\ r i =b + = i^| r=fe + since S' and S system hold the same coordinates in the domain 
r > b. Note that the boundary conditions Eqs. (jH]) certainly satisfy in virtual space. However 
the metric is discontinuous at the interface r = b in S system, therefore the relation of 
components Eg\ r=b + = Eg\ r=b - is not trivial, but a strong limit to the metric. 

By comparing the first and last term, we note that the electric field at r' = b is demanded 
to have a fixed relation to the field at x n '(b~, 9, (f>) in S' system. Because of the locality of 
electromagnetic fields, a certain relation of fields between different positions can not be 
established. Therefore (b, 9, 0) and x H (b~ , 9 , </>) must be the same point: 

r'(b-,9,4>) = b, 9'(b~,9,4>) = 9, 4>'(b-,9,cf>) = 0. (12) 

This is the requirement of non-reflection. Thus, we have 

(dff_ _ dt, 
\d9 ~ ~di 



dr' _ dr' _ 89' _ d<p' 
r =b- ~' \~d9 ~ Ikb ~ ~d6 ~ ~d9 



0. 



Substituting it into Eqs. (TTTj) . the equalities are directly satisfied. A similar study to the 
boundary conditions of magnetic field can give the identical result. 
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In fact, the interface at r = b in S system generally doesn't correspond to the interface 
r' = b in S' system for an arbitrary transformation. Thus the metric is usually discontinuous 
across the interface, and we even hardly measure the element of length on the interface. For 
example, if the inverse transformation is simply r' = f(r), 6' = 6, <fi' = as r < b, the 
elements of length on the inner and outer surface of the interface are respectively 

d/ 2 | r=6 _ = f(b) 2 (d6 2 + sin 2 £d0 2 ), dl\ =b+ = b 2 (d6 2 + 6 2 sin 2 M0 2 ). 

Nothing but f(b) = b gives d/ 2 | r= ;,- = d/ 2 | r= {,+ with arbitrary d9 and dtp. It is also available 
for cylindrical cloak by a similar analysis that non-reflection requires invariance of coordi- 
nates on the interface before and after transformation. 



III. ANALYTICAL EXPRESSIONS OF LIGHT-RAY AND WAVE-NORMAL RAY 
A. Spherical cloak 

Considering a spherical cloak, the cloak hides a homogenous-media sphere, and its center 
is at the origin of a spherical coordinate system. A plane wave E 1 ^ = E e lk ° rcos9 e x with 
x polarization is incident to the cloak. It has been proved that both the fields E 1 ^ inside 
the internal area (r < a), and the scattering fields E^ } outside cloak (r > b) varnish if the 
constitutive parameters of cloak satisfy the formula Eq. (j7]) and the boundary conditions 
f(a) = 0, f(b ) = b. Meanwhile the electromagnetic fields in cloak layer (a < r < b) have 



been given in 15] 

4 C P) = £o (f'(r) sinflcos0e r + ^cos#cos0e e - ^sin^ej e ifco/(r)cos9 , (13a) 
D c {p) = e E ^- (J^-sm9cos<f)e r + f(r) cosflcos0e e - f(r) sin0e^J e ifc( '^ cose , (13b) 

# c = ( j'( r ) sin # sin 0e r + li*l cosdsin ie + cos 0ej e ifco/(r)cose , (13c) 
rj \ r r J 

g = Eottr) ( f(d s i n s i n( f ) £ r + f( r ) c os9sm<f>ee + f'(r )c os0e A e ifco/(r)cos9 , (13d) 
1 ' cr \ r J 



where r] = y/io/^o an d k = w/c are the impedance and wave number of vacuum respec- 
tively. Since only fields in P space are discussed in this section, we omit the subscript "(P)" 
in all the expressions of coordinate, but still retain it in the expressions of electromagnetic 
quantities. According to above fields , we obtain the time-averaged Poynting vector 

(S° m ) = ^Re(E^xH^) = ^I^-fM cos 9e r - f'(r) sin 0e*Y (14) 
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— * 

t can be verified that the Poynting vector satisfies energy conservation V • (S? p) ) = 



23l |. Supposing the medium of cloak is non- dispersive, the time-averaged energy density of 



electromagnetic fields has the simple form 

IW< \ 1 kJ 1 F< n c fr\ e E 2 f(r) 2 f(r) 

Thus, the ray (energy) velocity in cloak is given as 23|, [24 1 

^ray = = C (j^ COS 6e r ~ sin O^j . (16) 



The ray vector is 24] 



s 



w \f(r) fir) 



In terms of the eikonal ip = k f(r) cosO, the geometric wave-front, ip = constant, is 

f(r) cos8 = const. (18) 
Then, we can get the wave vector fc (P) , which is orthogonal to the wave-front, in the cloak 

fc (p) = VV> = h ( f'(r)cos9e r - ^-sm6e e ) . (19) 



In the anisotropic media, the ray vector ( time-averaged Poynting vector ) is tangent to 
light-ray, however the wave vector is tangent to wave- normal ray, where s and k are generally 
not parallel to each other. The parametric equations of light-ray can be written as 

dr 1 

dA = a < p >" = kJV) C °* 9 ' (2 ° a) 
dA k f(r) 

^ = s m * = 0, (20c) 

where A is the parameter of light-ray parametric equation, and the components of ray vector 
employed above are contravariant in holonomic spherical coordinate bases. The <fi maintains 
its initial value 0o through out a light-ray in terms of Eq. (I20cl) . By eliminating A from Eq. 
pjaj) and Eq. fl20bl) , we have 



fir) J sin 6 

Then the analytical expression of light-ray is obtained 

f(r)sm6 = const, (j) = (f) . (21) 
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Here, the integral constant and 0n depend on the incident point on the outer surface of the 
cloak. 

A similar analysis is used to get the parametric equations for the wave-normal ray 

dr 



dA 
d9 

dA 

d0 

dA 



k 7 
k 6 
h < 



k f'(r) cos 6, 



An 



sin 9, 



0. 



Thus, the wave-normal ray equation is 



cxp 



-dr 



sin 9 



const, 



(22a) 
(22b) 
(22c) 

(23) 



r 2 f'(r) 

According to Eq. (I16p . when r closes to a, f(r) tends to 0, and the modulus of ray 
velocity increases beyond vacuum light speed c and tends to infinity which is against the 
law of causality. The contradiction indicates that the assumption of non-dispersive cloak 
is unpractical. In the presence of dispersion, the time-averaged energy density has a more 
complicated formula. It includes the derivative of constitutive parameters with respect to 
oj |24j, and then the super luminal result should disappear. More detailed discussion to this 
point can be found in 18|, ll9[. Nevertheless, the expressions in Eqs. ( IT6l [T71) still exactly 
represent the direction of the "real ray velocity" , therefore we will retain the expressions of 
V(p)ray an d s (P) , but limit them only to stand for the tangent vector of light-ray. 



B. Cylindrical cloak 

Considering a cylindrical cloak with indefinite length, its axis is z axis of a cylindrical 
coordinate system. The incident wave is a transverse-electric (TE) wave E 1 ^ = E e lk ° rcose e z 
with z polarization propagating towards x direction. If the cloak hold the perfect invisible 
parameters expressed in Eq. (jHJ), the calculation in Appendix B gives E 1 ^ = in the inner 
hidden area (r < a), and the scattering E^ } = outside the cloak (r > b). Meanwhile the 
fields in the cloak layer are 

4 C P) = E e ikof{r)cose e z , (24a) 
= £ ^o //(r) r /(r) e ifco/(r)cose e z , (24b) 

H° = -— (fir) sinfle r + ^ cos#e e ] e ifco ^ cose , (24c) 
Vo \ r J 

S ( c p) = sm9e r + f{r) cosOe^j e ikof{r)cosd . (24d) 
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Thus, the time-averaged Poynting vector is 

(4)) = ^(^cos9e r -f'(r)smee y (25) 

and the time-averaged energy density of electromagnetic fields under the assumption of 
non-dispersion is 

{Kr)) = S^Wr) (26) 

By calculation, the ray vector and wave vector in the cylindrical cloak have the identical 
form in the spherical cloak as shown in Eq. (TTT|) and Eq. (j!9p . Whereas, the variables r and 
9 in the spherical cloak are spherical coordinates, while in cylindrical cloak they represent 
cylindrical coordinates. 

Thus, following the same calculation as in spherical cloak, we derive the light-ray expres- 
sion in cylindrical cloak 

f(r)sm9 = const, z = z , (27) 

and the wave-normal ray equation 

" f f( r ) 1 

exp / dr sin^ = const, z = z . (28) 

[J r 2 f'(r) 

C. Examples of four types of cloak 

In this section, we apply the obtained expressions of light-ray and wave-normal ray to four 
cloaks constructed by different transformation functions. Type 1 is the linear transformation 
function fi(r) = [b/ {b — a)](r — a) provided in [l|, Q]. Type 2 is /^(O = b[(r — a)/(b — a)] 2 , 
and type 3 is /^(r) = b — b[(b — r)/(b — a)] 2 . Both type 2 and 3 are from [15]. Type 4 is 



/4(f) = (b/2a)[2a — b + Vb 2 — Aab + 4ar] provided in 16|. They all meet the invisibility 
condition f(a) = 0, f{b) = b (see Fig. 2). 

We construct four spherical cloaks with the four transformations. Each raw in Fig. 3 
displays a type of cloak corresponding to /i(r) to f±{r) from top to bottom. The first 
column displays the components of electric fields E (P) ^ in the x = plane. The second 
column displays the light-rays and wave-fronts of each cloak. The third column shows the 
wave-normal rays and wave-fronts which are always perpendicular to each other. 

From Fig. 3, some characteristics of each cloak can be concluded. In Fig. 3(b), the light- 
rays and wave-fronts hold equal spacing respectively, which is due to its linear transformation 
/i(r). For the second cloak, because f^ij) — > as r — > a, the 9 component of energy flow 
decreases in r — > a so fast that the energy traveling in cloak mainly concentrates at the outer 
boundary. Thereby the light-rays mostly compress near the outer boundary shown in Fig. 
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FIG. 2. Transformation functions. The inner and outer radiuses of cloak are a = lm and b = 2m 
respectively. 

3(e). For the third cloak, since f'^{b) = 0, the light-rays are condensed to inner boundary as 
shown in Fig. 3(h). In addition, r = b as k (P)r = 0, that is to say, the circle r = b becomes 
a closed wave-normal ray, and the wave-normal rays inside cloak are tangent to the outer 
interface r = b, as shown in Fig. 3(i). For the fourth cloak, its e and /z are continuous to 
vacuum across outer interface due to f'^b) = 1. Therefore, all the light-rays, wave-normal 
rays and wave-fronts are smooth when they pass through the outer boundary of cloak as 
shown in Fig. 4(k)(l). 



D. Verifying the light-ray equation is the solution of geodesic equations 



In this section, We will verify that the light-ray analytical expression (1211) in a spherical 
cloak satisfies the geodesic equation in S system of virtual space, which is [17[ 



^ + r$,^ = o, 



(29) 



where = lj im (j m j,i + 1mi,j — 1ji,m) is the Christoffel symbol in 3-D space. In S system, 
the spatial metric is 



7ij = diag(f(r) 2 , / 2 (r), / 2 (r)sin 2 #). 
Substituting it into Christoffel symbol, we have 



(30) 



rr 


f"(r) 
f'{r) ' 


-pr 


_ fir) -rr _ 
fir)' 1 4><t> 


fir) ' 




f'(r) 
f(r) ' 


r 9 


= — sin 6* cos 6* , 






fir) 
fir) ' 


p0 
1 cj>9 


= cot 9, 





(31) 
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FIG. 3. Lights travel in x = plane of four space where four types of cloak are placed respectively. 
Field distributions of E^i x \ of each cloak are shown in (a)(d)(g)(j). In (b)(e)(h)(k), blue lines 
denote light-rays, red lines denote wave- fronts. In (c)(f)(i)(l), green lines denote wave-normal rays, 
red lines still denote wave- fronts. 
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and other components are all null. Then, the geodesic equations become 

+ Q±kTk' - M*V _ M sin * k * k * = 0, (32a) 
dA f'(r) f'(r) f'(r) 

c\k e f'(r) 

^- + 2^4k e k r -sinO cos Ok+k* = 0, (32b) 
dA f(r) 

^+2^4^ + cotdk°k+ = 0. (32c) 
dA f(r) 

Tracing the ray path shown in Eq. (I2TI) . the differentials of coordinates satisfy 

f'(r) 

d9 = -^-fftan^dr, d<p = 0. 
/(0 

Substituting them into the relation between parameter A and element of length dl = 
(jij dx l dx j y/ 2 , which is dA = — {cy/— g 00 /u)dl [25], we have 

dA = fffrldr. 

CO COS 

Therefore, the wave vector can be expressed as 

dr _ cos# e _ d^ _ sin0 _ 
* - dA - fc o //(r) , k - ^ - fc 0/(r) , fc - U. (33) 

This is the wave vector in S system, its relation with the wave vector ([19]) in P space will be 
discussed in section 4.2. Substituting k % into Eqs. (1321) . then the left sides of the equations 
reduce to zero, equal to the right side. As a result, we find out that Eq. (|2ip is a particular 
analytical solution of the geodesic equations. 



E. Physical interpretation of light-ray equation 

We have verified the light-ray ( J2T|) denotes a geodesic expressed by curved coordinates S 
of virtual space, thus the physical meaning of the light-ray equation is quite clear. 

Since the virtual space is a flat vacuous space with a zero curvature, the light-ray trajec- 
tory which is also the geodesic in virtual space is just a straight line. Fig. 4(a) shows the 
light-ray with a Cartesian coordinate system (S' system). Through the inverse transforma- 
tion r' = f(r), where f(a) = 0, f(b) = b inside the domain r < b in S' system, we obtain a 
new system, S system. Through the transformation, the domain r' < a in S' system shrinks 
into origin, that is, the coordinates with component r < a do not exist in S system. Fig. 
4(b) shows the result in S system. The coordinate lines, whose tangent vectors are bases, 
are defined by 

Latitude lines : y = rsin^ = / -1 (r')sin#' = const. (34a) 
Longitude lines : x = rcos9 = /~ 1 (r')cos^ / = const. (34b) 
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FIG. 4. Interpretation of light-rays, (a) The red line denotes a light-ray propagating in virtual 
space. The meshes are the coordinate lines of S' system, (b) The red line is still the light-ray in 
virtual space. The curved meshes are the coordinate lines of S system. Two orange dashed lines 
are the boundaries inside which latitude lines shrink into the origin point, (c) Light-ray travels 
in P space with Cartesian coordinate meshes, (d) The mapping of lines with the expression of 
geodesic in S system to P space. 

where f~ l {r') is the inverse function of f(r). We observe that the coordinate lines cave 
towards the origin, especially the lines whose \x\ < a or \y\ < a outside the interface 
r = b cave directly into the origin when they pass through the interface. Therefore, the 
origin becomes an ambiguous point of coordinates in S system. However, the coordinate 
transformation doesn't change the curvature of virtual space, and the light trajectories, i.e. 
geodesies, are still straight with the expression in S system 

y = r'sinO' = f(r)sm9 = const. (35) 

It corresponds to the light-ray equation that we have obtained previously. 

The light-ray function maps to P space with the same form /(r (P) ) sin# (P) = const. Since 
P space is also a flat space, the light-ray becomes convex as shown in Fig. 4(c). Meanwhile, 
the constitutive parameters both in spherical and cylindrical cloak have singular components 
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with zero or infinite value at inner surface for f(a) = 0. In other words, the geometrical 
singularity of coordinates turns into the singularity of material, when mapping from S system 
to P space. 

We use an analogy to describe the process. The virtual space can be compared to a piece 
of flat elastic cloth with straight meshes on it. The meshes look like S' system. If we fix 
the boundary r = b, and drag the whole region r < a into the center point, the meshes on 
the cloth are curved as the coordinate lines of S system. Then, we draw a straight line on 
the cloth just as the red line in Fig 4(b) which plays the role of light-ray in S. After that, 
we loose the piece of cloth held into the center, and let it return flat. Then the meshes also 
returns to straightness, but the red line we have drawn becomes convex exactly as the ray 
propagating in P space. 

Now we should clarify the meaning of Fig. 4(d). Sometimes the meshes in this figure are 
interpreted as the coordinate of S system [4|. This interpretation gives rise to an illusion 
that S system is a real curved space, as well as the light trajectories. However, the meshes 
in Fig. 4(d) are neither coordinate lines nor geodesies no matter in S system or in P space. 
The only reasonable interpretation is that they are the curves in P space directly mapping 
from the straight lines in S system, and can describe the light trajectories in P space [3,0]. 



IV. DUALITY PRINCIPLE IN ANISOTROPIC MEDIA 

In virtual space, light propagates in straight lines, and the directions of wave vectors 
and ray vectors are always identical, however when mapping to physical space, they are not 
parallel to each other any more. So what exactly causes the split? Why does the light-ray 
instead of not the wave-normal ray in P space inherit its expression in S system? In this 
section, we will try to answer these questions by using duality principle. 

In charge free anisotropic media, electromagnetic fields have the following duality rule 



23] 



D E B H e jt k/u v p 

Duality rule 



E D H B f 1 /2 1 slo 1/u, 



ray 



where v p = u/k is the phase speed. In terms of the above systematic interchange relations, 



the either side satisfies a set of formulas for plane waves 24j: 
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D 
B 



e- E 



—k x H 
to 

-k x E 



to 



duality 



kjuj = k/v p 



Dx B 



E 
H 



jr l -b 



■ —LOS X B 

= usxD 
Ex H 



(36) 



E ■ D D-E 

According to the left side of the formulas, the wave equation is founded as 



k x 



/T 1 ■ (kxE) 



+ u 2 e-E 



0. 



(37) 



Concerning the special case of impedance matched material e/e = jl/ /i = n, the dispersion 
relation, which are also called wave- vector eikonal equation, is (derived in detail in Appendix 

c)B 

(38) 



k i k j - 1 {^j det(n^) 



2 1 fu\ 2 

nijk l k J ( — ) det (rii 



0. 



7 V c / 

Here, there are two ambiguities to be clarify. First, the expression given in [4] is based on 
Cartesian coordinates, so 7 reduces to unit, while the expression here is valid in any curved 
coordinates. Second, n u and presented above express the components in holonomic 
coordinate bases (Landau's definition), however, the dispersion relations expressed either in 
Minkowski's definition or in anholonomic unit basis are a little different from Eq. ( |38l) (see 
Appendix C). In [4], the three kinds of expression are congruent with each other in Cartesian 
coordinates. 

Using the duality principle, we can easily obtain the dual form of wave equation 



s x 



jl- (sx D ) 



D 



0. 



Thus, the "ray-vector eikonal equation" is 



in 



(n-%s % s J ~ - {-) det l(n-%] 



0. 



(39) 



(40) 



where (n is the inverse matrix of n^, and it is easy to verify that its corresponding 
covariant components {n' 1 )^ are also the inverse matrix of n* J . 



A. Eikonal equation in curved coordinate system S 

We first consider the relation between k % and s % in curved coordinate system S. In S 
system, the relative constitutive parameter is n u = (n" 1 )*- 7 = 7^. Substituting it into Eq. 
( 1381) and Eq. ( f40|) . we get the wave- vector eikonal equation in S system 

u\ 2 , ■ • /u\ 2 



^hk 3 - (-) = % .AW-(-j =0, (41) 
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and the ray- vector eikonal equation 

j'siSj - (£f = - = 0. (42) 

The two equations indicate that \k\ = k and |s| = l/k . By considering the general 
relation k ■ s — ^iijk % s^ = 1 [24] . we can conclude that the directions of k and s completely 
coincide, and they have the relation fcj = fco 2 -^- This result is actually not marvelous, since 
the electromagnetic wave is merely propagating in a vacuous space, despite the space is 
expressed in a curvilinear coordinate system with non-unit metric 7^. 



B. Eikonal equation in physical space P 

In physical space, the media of the cloak are inhomogeneous and waves traveling in it are 
no longer plane waves. Nevertheless, we can verify the fields obtained both in the spherical 
and cylindrical cloak still satisfy Eqs. (JSEl E3 EH])- The relative constitutive parameter of 
the cloak is n (P) = 1^/7/7^) 7 iJ , and its inverse is (fiT P \)ij = \/7(p)/7 Jij- Substituting them 
into Eq. (1381) and Eq. (|40l) . we obtain the wave- vector eikonal equation in P space 

7%)j*(p)/-(^) = (T(p)ife7(p)i/T W ) WW - (^) = 0. ( 43 ) 
and the ray- vector eikonal equation in P space 

7«%)V-(^) = (7(p)VSi)V)^(p) J -(9 = °- ( 44 ) 

Comparison of Eqs. (1431 144)) and Eqs. (T4Tj |4*2|) manifests that the wave-vector eikonal 
equations expressed by covariant components ki have the same form in P space and S system. 
However if we write the wave- vector eikonal equations with contravariant components k % , 
the expressions in P and in S are different. Similarly, the expressions with contravariant 
components s l have the same form in the two systems, yet the expressions with covariant 
components Si are different. Considering the relation k (P) ■ s* (P) = 1 simultaneously, we can 
get 

k ( P H = h 2 7yS (P) l , S(p/ = k ~ 2 7 lJ k (PH . (45) 

This relation is easily validated according to Eq. ifPTj) and (fl9l) . It is also valid in S system 
which means k is parallel to s, however this does not make sense in invisibility cloak. 

These results also can be verified with canonical equations given by [4]. If we define 
Hamiltonian by wave- vector eikonal: 

H = ^k (P)i k (P)j -(^)\ (46) 
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and according to the first equation of canonical equations 



<H>) = Mo* = (47) 

dr dk {P)i ' dr dx l {p) ' 

where r is the parameter of canonical equations, we obtain the tangent vector of light-ray 
dx* p) /dr = 2yik (P) j. It reveals the tangent of light-ray is exactly towards the direction of 
s (P) l given in Eq. (|45l) . If we define another Hamiltonian by ray- vector eikonal: 

H = (WW^Kr^CP),- - (£) , (48) 
and modify the canonical equations into 



dx* P) OH ds (P)i dH 



dr dS( F) i dr 9xl p) 



(49) 



we would have 

dx l 

= 2 7(P / fc 7 (P /W (P ), = 2A: - 2 ^ P) . (50) 

As a result, the modified canonical equations depict the trace of wave- normal rays. 

Actually, if we express ki and s 1 with D'\ B l , Ei, Hi in Minkowski's definition, we have 

, e ljk DiB k . le^ k D 3 B k 

ki = to ' — — , s = — -=- — . (51) 

The expressions apply to both S system and P space. According to the relations in Table 
[H we can get two new corresponding relations between S and P: 



(p)- 



We note that k % and k l , p , do not correspond to each other, and neither do Sj and s (P) j. In 
other words, covariant components of wave vector k{ in S system map to wave vector k (P) i in 
P space, yet contravariant components of wave vector k l (= k^s 1 ) map to ray vector s l (p) in 
P space. Although k i} s l denote a same orientation in S system, k (P)i , and s l (p) denote two 
different orientations in P space. In addition, the contravariant components s l = dx l /d\ 
directly represent the directional change of coordinates dx l , therefore light-ray equation in 
P space inherits itself's expression in S. By contrast, ki can not directly evince the change 
of dx*, so the expression of wave- normal ray takes a new form in P space. 

According to the dispersion relation (|43p . we can obtain the group velocity in the cloak 
through applying the formula v g = Vgw, where is the divergence operator in fc-space. 
Supposing the medium is non-dispersive, the group velocity takes the identical expression of 
ray velocity ( fi~6i) . It means the group velocity and ray velocity are congruent in the cloak. 

As commented in section 3.1, the contradiction of super velocity of light takes into account 
the existence of dispersion. However, the dispersive processes of e and \i are different, so 
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the case of impedance matching established in a particular frequency would hardly retain 
when frequency changes. Therefore, the dispersion relation of impedance matched material 
(J55|) is unavailable when we use v g = Vgw to calculate the group velocity in dispersive case, 
because Eq. (|38|) can not display the differences between the derivative of e and fi with 
respect to u. 

V. DESIGN OF THE INVISIBILITY CLOAK WITH RAYS SATISFYING HAR- 
MONIC FUNCTION 

In Ref. |2j, Leonhardt applied }f\yKOBCKHH function = z + 1/z, where z = r e ie , 
to a conformal mapping, and obtained an inhomogeneous but isotropic medium which can 
realize an invisibility device in geometric limit. In the media, the families of light-rays and 
of wave- fronts are the imaginary and real part of ((z) respectively. In addition, they are 
both harmonic functions and orthogonal to each other. 

Using the expression of light-ray (1271) . we can design an invisibility cloak by means of 
transformation-optics method, which gives the same light trajectories as in Leonhardt's 
device. However, we start our work with a more general target to find the general form 
of transformation function f(r) which causes the light-rays in cylindrical cloak to satisfy 
harmonic function. Substitution of light-ray (J27j) into Laplace's equation yields 



Then, it reduces to Euler's differential equation, r 2 f"(r) + rf'(r) — f(r) = 0, therefore the 
general solution is f(r) = C\[r 2 — C^jr . The function has a zero point r = C*2 1/2 > 0, which 
is the most important condition to construct an invisibility cloak. If we set C\ = b 2 / (b 2 — a 2 ) 
and C*2 = a 2 , the transformation function becomes 



The result evinces that each solution corresponds to an invisibility cloak with a spacial 
pair of inner and outer radii. In the process of derivation, we didn't command f(a) = 
0,f(b) = b, however the result indicates that the form of harmonic functions V(r, 6) = 
f(r) sin 6 is unique, and it always can realize invisibility cloak. For convenience, we omit 
the coefficient of f(r) in the following derivation which doesn't influence the shape of rays. 

V = const represents the family of light-rays in the cloak. Let V be the imaginary part 
of an analytical function F(z). Using Cauchy-Riemann condition, we obtain the real part 
U(r,6) = \{r 2 + a 2 )/r] cos 9. The U = const stands for a family of curved surfaces, which is 
orthogonal to the energy flows. The analytical function F(z) is 



V 2 _ D [/(r)sin0] 



0. 



(52) 



F(z) 




+ — . 

z 



(54) 
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FIG. 5. Cloak with rays satisfying harmonic equation, (a) Blue lines denote light-rays, orange 
lines denote the surfaces orthogonal to light-rays. In cloak layer, they are painted by Eq. (|54p . (b) 
Green lines denote wave-normal rays, red lines denote wave-fronts. In cloak layer, they are painted 
by Eq. flggp . 

If we simply set z' = z/a, then F(z) = a((z'), where ((z) = z + 1/z is precisely /Kvkobckhh 
function. It is exactly the conformal mapping which Leonhardt used to design the invisibility 
devices 0]. The complex analytical function is a compact way to describe the ray path by 
the meaning of real and imaginary part explained above, as is shown in Fig. 5(a). 

Regarding to the transformation (!53|) . the eikonal ip = kof{r) cos 8 is a harmonic func- 
tion as well. So we also can compress the wave-fronts and wave-normal rays into another 
analytical function 

G(z) = z - — . (55) 

The real part represents the family of wave fronts, and the imaginary part represents the 
family of wave-normal rays exhibited in Fig. 5(b). 

Owing to the energy conservation, the time-averaged Poynting vector (S? p .) is divergence- 
free. And the wave vector must be irrotational, since /c (P) = Vip. However, if and only if the 
transformation satisfies Eq. ( 1531) . we can prove both (S? p) ) and k (P) are not only divergence- 
free but also irrotational. This conclusion is consistent to the fact that both light-rays and 
wave-normal rays satisfy Laplace's equation for the particular transformation function. 



VI. CONCLUSION 

In conclusion, we have meticulously compared the behaviors of electromagnetic fields in 
P space with in S system of virtual space, and provided the corresponding relations with 
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electromagnetic quantities and boundary conditions between the two systems. Then, we 
obtain the analytical expression of light-rays and wave-normal rays, and give their physical 
interpretation. Furthermore, we verify the expression of light-rays we obtained is exactly a 
particular solution of geodesic equation in S system. In addition, using duality principle, we 
show that the covariant and contravariant components of wave vector in S system map to 
wave vector and ray vector respectively in P space, and it causes the split of k (P)i and s* p) in 
P space. Finally, we find out the general form of transformation function which turns the 

— * — * 

light-ray to harmonic, and therefore we can construct a special cloak in which (S^ p) ) and k (P) 
are both divergence-free and curl- free. 
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APPENDIX A: DIFFERENT WAY TO DEFINE ELECTROMAGNETIC FIELDS 
IN CURVED SPACETIME 

In Minkowski spacetime, 3-D macroscopic Maxwell's equations in a Lorentz frame are 

diB' = 0, ^- + e ^ k d E k = 0, (A.la) 

diD* = p, -— + e i i k d j H k = 0, (A.lb) 
and their 4-D covariant forms are 

dxF^ + dvF^ + dpFrt = 0, (A.2a) 

d a G aP = J 13 . (A.2b) 
The four 3-D electric and magnetic vectors, E iy H iy D l and B % are defined as 

Ei = F 0i , B l = ~^ k F jk , D l = ~G 0i , ^ = ~e ijk G jk . (A.3) 

Here we use Greek indices a, (3, - ■ ■ to run from to 3 for time and spatial parts, and stipulate 
the sign of metric takes the form ( — , +, +, + ). For vacuum, G a(3 = cE F al3 . 



23 



In curved space time, 4-D microscopic Maxwell equations which have the same form as 
the macroscopic Maxwell's equations in vacuum become 2l|, |25] 



dxF^ + d v F\fi + d»F u 



ce Q V a F°* 



CSq 



-gF 



a(3\ 



(A.4a) 
(A.4b) 



where V Q is 4-D covariant derivative, g is the determinant of the spacetime metric g^ v . 
If the excitation tensor G Q/3 is defined in the form 0, Q, S, 0] 



G aP = ceoV^gF 



a/3 



(A.5) 



and the definitions of 3-D electric and magnetic vectors are retained for flat Minkowski 
spacetime as in Eqs. (IA.3j) . both 3-D and 4-D Maxwell's equations keep the form of flat 
spacetime in Eqs. 1 1A.1I) and Eqs. ( 1A.2I) . except the 4-current charge density p, and 
3-current j l should be replaced with the reduced forms J 13 = \f—gJ^, p = a/7P, f = \fl? 
respectively |7|, where 7 denotes the determinate of spacial metric 7^ = g^ — 9oigojl9oo- 
We name this set of excitation tensor and 3-D electric and magnetic vectors as Minkowski's 
definition. 

The constitutive relations in Minkowski's definition are 



D l 



B' 



e 



Po 



7 



"000 



1 



Jjk 



cgoo 



gojHk 



7 



-900 



Jjk 



cgoo 



gojEk-, 



(A.6a) 



(A.6b) 



Thus, the permittivity and permeability satisfy = £o\/—l/gool l:i , P %i = Po\/— ll gool l K 
and magnetoelectric coupling tensor is = e t:jk g j/(cg Q). If g^ v satisfies g 00 = — l,<?oi = 0, 
then magnetoelectric coupling vanishes, while and p l i reduce to Eqs. (@J. 



Another definition of 3-D electric and magnetic vectors appears in Landau's book 25], 
therefore, we call this set as Landau's definitions. In Landau's definition, we put a bar on 
each quantity, and the excitation tensor retains its formula in flat spacetime G al3 = ceqF 
and the four electric and magnetic vectors are redefined as 

1 



25 



28] 



E; 



0i> 



B' 



2 Cv /7~ ~ jkl c ' ' 2 

Under this set of definitions, 4-D Maxwell's equations have the form of 

d x F^ + d u F Xfl + d^x = 0, 

\7 n G a(S = J 13 . 



D' 



y-9oo QQj 



Hi 



-9 



e ijk & k . (A.7) 



(A.8a) 
(A.8b) 
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The corresponding 3-D Maxwell's equations are 

-^(Vt^) = 0, -LJL(^S') + ^%E k = 0, (A.9a) 

1 d :(y/jD i ) = p, -J-O^D*) + ^%H k = X. (A.9b) 



Comparing Eqs. (1A.3I) with Eqs. (1A.7jl . we have G alS = y/—gG al3 , and the transformations 
of other quantities between the two definitions are identical with Eqs. (J2j) . 

Assume that an observer whose 4-velocity is u a in a curved spacetime endowed with an 
electromagnetic field F^ u . Then, questions arouse naturally. What is the measurement of 
electric and magnetic fields by the observer? Is the measurement equal to one of the two 
definitions mentioned above? Actually, there has been a mature way to answer this question. 
In the curve d sp acetime, electric and magnetic field vectors are defined respected with the 
frame u a as [29] 

E a = -F aP u p , H a = *G a pu 13 = X -e a?llv G^vP, (A.10) 

where *G a p = € a p flu G tJ ' u /2 is the dual of G^ v . Here, we put a tilde on electric and magnetic 
vectors to discriminate its representation under Minkowski's definition and Landau's defini- 
tion, however, the definition of excitation tensor retains its formula from Landau's definition 
Qa/3 _ c£fj p&l3 ^ w hich guarantees that G is a tensor of 4-dimension. 

In terms of the definitions of Eqs. ( lA.lOj) . E a u a = H a u a = 0, which evinces that E 
and H are both space- like vectors for the observer u a . In addition, if we simply let electric 
displacement D a = EoE a and magnetic induction B a = fioH a , i.e. what are made in vacuum 
of flat spacetime, the components of E, H, D, and B in the tetrad e ( ° ) of the observer, where 

e (0> = u0l y are 

= F {0l) , H {1) = -~e ijk &i k \ m = -- c G {0l) , S<*> = --e^F^y (A.ll) 

As we see, the components in tetrad satisfy the original relation of flat spacetime (1A.3I) . 
Furthermore, the electromagnetic tensors admit the following decomposition 29^ 31 1 : 

F^u = (E^u u - E v u^-ce^ aP u a B^ (A.12a) 

giv = c {D^u u - D u u^) - e^UaHp, (A. 12b) 



and a further study reveals the form of Maxwell's equations expressed with E and H [30_|, |31|. 

Here, we are interested in the condition that the observer is comoving. For a comoving 
observer, u a = (l/^—g 00 , 0, 0, 0), and the transverse projector h^ u = g^ u + u^u v reduces to 
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the spatial metric 7^ mentioned earlier. In this case, the relations ( 1A.12|) lead to the result 



G jk 



-c^fye ijk B 



V-900 



7^ E ^ 9k ^ 



(A.13) 



It seems that the measurement of electric and magnetic fields by an observer conforms 
neither to Minkowski's definition nor to Landau's definition. However, if we set goo — 
— l,<7oi = 0, the reduced relations Eqs. (1A.13I) are identical with the reduced relations Eqs. 
flA.7j) under Landau's definition. So, if we use artificial materials to simulate the behavior 
of electromagnetic fields in curved spacetime based on the form equivalence of Maxwell's 
equations, it should be heeded that the distinction between the electric and magnetic fields 
measured in materials and measured by the observer moving in real curved spacetime. 



APPENDIX B: ELECTROMAGNETIC FIELDS IN CYLINDRICAL CLOAK 

In this portion, we would only consider TE polarized wave. Substituting the constitutive 
parameters Eq. ©into wave equation fl37|) . we get the general scalar wave equation for 
all transformation /(r) which is exactly the Helmholtz equation in original S' coordinate 
system of virtual space: 

1 d ( dE (P) \ 1 d 2 E (P) 9 , 

mV^f) + p^ + kE "" = - (B1) 

Through separating variables, electromagnetic fields can be represented by series 

00 

E (P) = Yl l a nJn{%) + KN n (x)] e in9 e z , (B.2) 
—00 

where J n (x) and N n (x) are n-order Bessel function and n-order Neumann function respec- 
tively. These expressions are available in every area, however x takes different forms in each 
area. Outside the cloak (r > b), x — k r is in the expressions of incident wave E 1 ^ and scat- 
tering wave E^p y In the cloak layer (a < r < b), E? p , takes x = kof(r). Inside the internal 
hidden area (r < a), x — k x r is in the expression of E™), where k\ = u^/eifii and £1, Hi are 
the material properties of the hidden medium. The incident wave E^ } = E Q e lk ° rcose e z can 
be expanded as 

00 

= EoY^ i" J n (k r)e ine . (B.3) 
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Because E sc should tend to an outgoing wave as r approaches to infinity, it is more convenient 
to expand scattering wave with Hankel function of first kind: 

oo 

4 S p, = E «Me W . (B.4) 

— oo 

In addition, = means that fields should be finite when r — > 0. Then applying the 
continuity of E z and Hq at inner and outer surface of cloak, we have following equations: 

a^Uha) = a c n J n (k f(a)) + b c n N n (k f(a)) , (B.5a) 

at-hUha) = ^[<j;(WW) + b c n N' n (k f(a))}, (B.5b) 

E i n J n (k b) + a%HV{kab) = < J„ (**/(&)) + b c n N n (k f{b)) , (B.5c) 

E Q rJ> n {k Q b) + < C #«'(M = M[a c n J' n (k f(b)) - b c n N' n (k f(b))} , (B.5d) 
Solving the group of equations, we obtain 



(B.6a) 



< = £ i" [J n (k b)HV'(k b) - ^(* to 6)irW(fc b 6)]/{ir«'(fc b 6)[J B (*b/(6)) 
+A^ n (W(&))] - MHVfab) [J' n (k f(b)) + A n N' n (k f(b))] 



aL nt 



/foofciJ ra (fc /(o))J4(fcio) - fJ.ikof(a)J^(k Q f(a)) J n (fcia) fl c _ ^ fl c ( B gb ) 
^ 1 k f(a)N' n {k f{a))j n {k 1 a)-^ak 1 N n {k Q f{a))j' n {kia) 

-rjT — r + A, N n (k f(a)) } <, (B.6c) 

< = ^)L {^[ J ^ W(&)) ~~ ^K(WW)]< -£ i n j;(M)}- (B.6d) 

Substituting invisibility condition /(a) = 0, /(&) = 6, we have = i?o i n , = ( i.e. A n = 
), therefore a^ c = which indicates no scattering. However, when calculating a™\ there is 
a term A n -N n (k f(a)), in which A n — > while N n (k f(a)) — >■ oo as /(a) — )■ 0. A meticulous 
calculation gives the result a™ 1 — > 0, so waves can not spread into the hidden area. Finally, 
we get the expressions of electromagnetic fields in the whole space as shown in Eqs.( l24l) . 

A detailed derivation should be identified. If we substitute the invisibility condition 
f(a) = 0, f(b) — b at the beginning, the fields in the cloak layer must have no terms of 
Neumann function, which means b c n = , since N n (kof(a)) — > oo is against the requirement 
that fields should be finite. Therefore, the last terms vanish in Eqs. (1B.5aj) and ( 1B.5bj) . In 
this case, the Eq. flB.5bj) still leads to a™* = 0, yet the result can not satisfy Eq. (IB.5aj) 



when n = 0. In the more general derivation above, when n = 0, though b c n — > 0, the product 
b c n ■ N n (kof(a)) is towards to a finite quantity —E i n /J n (kia), and it counterbalances the 
equality flB.5aj) . however, the product only exists at the inner surface, so it no longer has 
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;he meaning of fields but represents the magnetic surface current flowing along the z axis 



3. 



APPENDIX C: EIKONAL EQUATION OF DIFFERENT FORMS 



Wave equation ( 13 7p can be written as the components form 



riii 



The condition which protects E l from trivial solution is 



E l = 0. 



det 



Expanding the determinant 



0. 



det 



e ipj k p {n x y k e kq ik q + ( — J n a 



on - 



3 r 



n 



a=l 



e i a p a j a k Pa (, n i y C " "tkaqaak 9 ™ + ( — ) Ui a 



UJ 



*2*3 



,a=l 
3 



+ (-) Yl Mi 

a C 



a=l 



a=l 



(C.l) 



(C.2) 



(C.3) 



where Mjz = e ip jNf^ = e ip jk p (n 1 ) J ' e kq \k q ^ and the indices of e ia j aka satisfy a mod 3. 
Calculating each term in the expansion, we have 



1st term: \{ M laOL = y/y^ei^e^e^N^N^N^ 



a=l 

f — f ■ f ■) N pijl N P2j2 N P3j3 — n 

PlP2j2 C JlP3j3 t JlP2j2 t PlP3j3/ iV l JV 2 JV 3 ~~ U ' 



(C.4) 



2nd term: 



(7)" E t in2l3 ^M la{a+l) M la(a+2) 



a=l 
3 



(7) E ^ la+lla+2 n laa M la {a+1) M la (l , 2 , 



a=l 
3 

£ 

a=l 



OL £p a + 2 ja+lja+2 



KJ-Pa + lja + l A j-p a + 2ja + 2 
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because 



2nd term • det (n^) 



= V / 7(^) 

a=l 

■ "A^(n- 1 ) i ^^c^ fa+a( a +2 )* tota 
= l(^y[^i a ant p pk ia k t Pjk a k fi ] = (j^fiijk^ 
therefore we have 



[^7e* 1 * 2 * 3 n tl in t22 ni3 3 ] 



2nd term = ( 7— riijtfk 3 J / det (ny) 



(C.5) 



4 

3rd term: (7) E e ™ M ^^( 



(a+l)"-i Q+2 (a+2) 

a=l 
4 3 

/7 u) E' 
0=1 

w 4 3 

(-) Yl [^(«+i) e 9 Q a(a+2) - V(«+2)e to a(a+l)] k Pa k Qo 
C a=l 



2- (^y 1 ,/^///,-'. 



(C.6) 



4th term: 



n 



7?.; 



det (n. 



(C.7) 



Substituting them into Eq. (lC.3p . we have 

det e ipj k p {n 1 ) jk e kql k' 1 + (^f n a = ^ff.) - ^ (^'det (%) = 0. 

7 (C.8) 

This is exactly the wave- vector eikonal equation ( l38l) expressed by Landau's definition. The 
ray-vector eikonal equation (1401) can be simply obtained through duality rules. In terms of 
Eq. (|2bl) and Eq. (jHJ), we can get the two kinds of eikonal equation written by Minkowski's 
components and anholonomic components. 

Under Minkowski's definition, wave-vector eikonal equation is 

2 



n ij kikj - f-Vdet (//'') = //,,/// - 



c / 



and ray-vector eikonal equation is 



det(ny, 



0. 



(C.9) 



1 2 



det [(n- 1 )^] 



00 



det [(rr%-] = 0. (CIO) 
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The wave-vector eikonal equation expressed by the components in anholonomic unit basis is 
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